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Abstract 

Cosmic rays are the highest energy particles available for our study and as such serve as excellent 
probes of the effects of Lorentz Invariance Violations, which are expected to increase with energy. 
This general paradigm is investigated in this paper by studying the effects of such violations within 
the Coleman-Glashow model in which each particle species may have its own maximum attainable 
velocity even exceeding that of light in vacuo. The particular focus here is that the muon neutrino 
may have the maximum speed exceeding that of light. We show that such an assumption leads 
to the elongation of the decay lifetime of the pion that increases with energy over and above the 
time dilation effects. We provide a transparent analytical derivation of the spectral intensities of 
muon neutrinos and muons generated in the Earth's atmosphere by cosmic rays. In this derivation 
we not only account for elongation of the pion lifetime, but also for the loss of energy by the 
neutrinos by radiation of the electron-positron pairs through the Cohen-Glashow process, during 
their propagation. We then compare the theoretical spectra with observations of neutrinos and 
muons from large instruments like IceCube and BUST to set a limit of ~ 1CT 13 on the fractional 
excess speed of neutrinos over that of light. We also show that the ratio of the spectral intensities 
of downward and upward moving neutrinos at various angles constitute a diagnostic exclusively for 
the Cohen-Glashow process, which may be searched for in the IceCube data set. We conclude the 
paper with several comments, including those related to improvements of these tests when definite 
signals of GZK neutrinos will be observed. 



I. INTRODUCTION 



The study of several exciting aspects of high energy astrophysics and indeed of many sub- 
tle aspects of basic physics have been given a boost by the commissioning of large detec- 
tors of cosmic ray secondaries such as ANITA, IceCube, Auger, BUST, Kolar Gold Fields, 
Kamiokande and other experiments with collecting powers of ~ 100 km 3 [THT5]. These de- 
tectors have already detected ~ 10 9 cosmic-ray muons of median energy ~ 2 x 10 4 GeV and 
~ 10 4 neutrinos that allow the spectra to be determined up to ~ 10 6 GeV. The physics input 
regarding high energy nuclear interactions from accelerators, colliders and other sources help 
in reliably modeling the propagation of cosmic-rays through the atmosphere and qualitatively 
account for the observed spectral intensities of the muons and the neutrinos. Comparison 
of these spectral intensities with the model predictions then allow us to probe into primary 
cosmic ray composition at high energies and search for effects due to new physics, such as 
small violations of the Lorentz Invariance that may manifest themselves only at the highest 
energies. This paper is devoted to such an exercise. 

Violation of Lorentz Invariance is studied from two distinct perspectives. The first is 
exemplified by Michaelson-Morley and Hughes-Drever experiments which test the existence 
of preferred frames of reference and the anisotropy of Machian type long-range interactions 
of matter in the laboratory with astronomically distant matter. The extraordinary accu- 
racy achieved in such interactions validated relativistic theories of gravity, especially GR 
[T6t [T7] . The other perspective is exemplified by the theoretical considerations of Coleman 
and Glashow [TBI EE] > who accept the possible existence of a preferred frame, such as the 
frame in which the dipole anisotropy of the universal microwave background at 2.7K vanishes. 
In the preferred frame, the laws of physics are assumed to be invariant under translations 
and rotations. However they investigate the possibility that different particles could have 
maximum attainable speeds different from that of light in vacuo, and these speeds could, 
in principle, exceed that of light by a small amount. Coleman and Glashow have developed 
a perturbative framework to discuss the violations of Lorentz Invariance (LIV) with terms 
that are renormalizeable and are invariant under the SU(3) x SU(2) x U(l) gauge symmetry 
of the standard model. Going beyond the standard model, Kostelecky and collaborators 
have carried out extensive analysis of models where Planck scale physics yields signals in 
the propagation of photons, neutrinos and other particles that have the potential for being 
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observed in present-day or future experiments [201-22J. These later papers provide a compre- 
hensive overview of the physics and the observational status of these models: Observations 
of high-energy gamma rays from distant astronomical sources have also been used to set 
lower bounds on the energy scale at which quantum gravity effects lead to increase in the 
velocity of light with energy [23]. The aim of this paper is to discuss the bounds on LIV 
derived from cosmic-ray observations based on the formalism developed by Coleman and 
Glashow [T8| [19]. In this context, we may refer to the elegant review of earlier work by 
Bietenholz [21]. 

We begin by recalling briefly the earlier efforts in the field of cosmic rays to search for 
the effects of superluminal velocities. An excellent review of the efforts to observe tachyons 
[25] in cosmic ray showers is provided by R. W. Clay [26] . The air-shower group of the 
Tata Institute of Fundamental Research pioneered these studies by searching for energetic 
particles that arrive at the air-shower array some 10-50 fis before the main shower front of 
electron-positron pairs, muons and gamma rays initiated by cosmic ray particles of > 10 6 
GeV [271128]. 

Following the lead given by Coleman and Glashow [18J, with specific reference to the 
present paper, the early bounds on LIV using horizontal air showers were obtained by 
Cowsik and Sreekantan [2S]; detailed comments on this paper may be found in the papers 
of Coleman and Glashow [HJ EE] and of Halperin and Kim [30] . This later paper maps the 
violations of Lorentz invariance into violations of the Equivalence Principle. In a subsequent 
paper, Cowsik et. al. [31] have investigated the possibility that if similar effects can induce 

—> v e + l, then such a rate is far more strictly bounded. Stecker and Glashow [32] discuss 
the bounds on LIV of electrons based on observations of energetic cosmic rays. Similarly 
Stecker and Scully [3"3"l I3"4] have put bounds on LIV in the hadronic sector by consideration of 
the GZK cut-off [351 136]. Direct observations of the neutrinos from supernova 1987A [371138] 
allowed Stodolsky [39] and Longo [40J to set bounds on any excess speed of neutrinos over 
that of light at the ~ 10~ 8 level, a significant improvement over early results at accelerators 
[H]. In the context of OPERA experiments [411450] several ideas of interest have been put 
forward and we reference a few for completeness [2H1 I2S1 ESI El I5TH54"] . 

The particular focus here is to provide an analytical calculation of the spectral intensity 
of muons and muon neutrinos arising from the decay of pions produced by cosmic rays in the 
Earth's atmosphere. In carrying out these calculations, we have included the enhancement of 
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pion life-time and the decrease in the average energy transferred to the neutrino in pion decay 
due to any posited superluminal motion of the muon neutrino. Secondly, we have included 
the effect of such a neutrino losing energy by emitting electron-positron pairs during its flight 
even through vacuum, as pointed out recently by Cohen and Glashow [65] in the context of 
OPERA experiments. 

We will not embark here on the ambitious program of getting the best limits for the rich 
variety of the LIV modifications for the various particles involved. Instead what we will 
attempt here is more limited and yet clearly illustrates the potential reach of this approach. 
Here we will focus on the effects of modifying the energy-momentum relation for only to 
E v = p u (l + a), as suggested by Coleman and Glashow for modeling violations of Lorentz 
Invariance. Our analysis presented here exclude values of a values down to ~ 10~ 13 . This 
is achieved by providing a transparent analytical calculation for the propagation of cosmic 
rays in the Earth's atmosphere that accurately reproduces the known data when no anomaly 
is assumed, and then comparing the theoretical spectra for various values of a with the 
observational data. 



II. CALCULATION OF THE SPECTRAL INTENSITIES OF NEUTRINOS AND 
MUONS GENERATED BY COSMIC RAYS IN THE ATMOSPHERE. 

The earliest calculations of the fluxes of neutrinos and muons in the earth's atmosphere 
were due to Volkova and Zatsepin in 1961 and Zatsepin and Kuzmin in 1962 [661 EZ]- This 
was followed by a slightly more detailed calculation by Cowsik and collaborators in 1963 and 
in 1966 [68l 169]. This later paper also describes the experimental aspects of the detection of 
these energetic neutrinos with detectors located deep underground. Since then, the calcula- 
tions have progressively improved with the explicit inclusion of the inelastic cross sections 
for the production of pions and other particles measured with particle beams at accelerators 
[TUHT4"] . In this section, we derive analytical formulae for the spectral intensities of muons 
and neutrinos arising from the decay of pions and include the effects of posited superluminal 
speeds for the muon neutrinos. There are basically two effects: (1) a progressive lengthening 
of the pion lifetime [75J and the reduction in energy transferred to the neutrino in the decay 
process due to LIV effects [76] and (2) the loss of energy suffered by the neutrino during 
propagation owing to the emission of electron-positron pairs through the Cohen-Glashow 
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process |65j. We will describe these two effects below and present an analytical calculation 
of the cosmic ray fluxes. 

A. Kinematics of pion decay with superluminal neutrinos. 

In the discussion of the kinematics of pion decay we make the minimal assumption that only 
the muon neutrino has a maximum attainable speed exceeding that of light in vacuo, and 
work within the framework of the Coleman-Glashow model for LIV [TSJ [T5]. The analysis 
below follows closely our earlier work in relation to OPERA results [75] . The key assumptions 
for this analysis are the following: (1) The relation dE/dp = v, the velocity of the particle, 
(2) energy-momentum conservation holds, and (3) the positivity of energy for free particles, 
which excludes tachyons. Since the mass of the muon neutrino is in the sub-eV domain and 
our considerations are limited to neutrinos generated by cosmic rays at high energies, say 
above 10 GeV, we may safely neglect the neutrino mass and write 

E v = Pu (l + a) (1) 

where a is the superluminal parameter, a very small quantity with a C 1. Note that a = 25, 
where 5 is the LIV parameter defined similarly by Coleman and Glashow [T%| [TU] . 

The superluminal energy-momentum relation in Eq. (flj) suppresses the pion decay both 
through its effect on the matrix element of the decay and through kinematic effects, which 
become progressively more severe with the increasing energy of the pion. We begin with the 
description of the kinematic effects: The pions and muons follow the standard mass-energy 
relation 

£ i = (p2+ m 2)i ( 2 ) 

It is convenient to express the momentum four-vector of the particles as 

p n = (E n ,p n , 0,0), pf, = (E^p^up^O) and ^ = (E u ,p v i,p ut ,0) (3) 

where the subscripts I and t refer to the longitudinal and transverse components. We 
explicitly satisfy momentum conservation by choosing 

Vvi = VPtv, Vid = (1 - v)Ptt and p vt = -p^ = p t (4) 

The equation for conservation of energy now reads 

{pi + ml)^= [pl{\ - V ) 2 + p 2 t + + [pl V 2 + pi] 1 (1 + a) (5) 
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At cosmic ray energies all the momenta are large compared with the masses of the particles 
and the square roots in Eq. ^ may be expanded keeping only the leading terms. This leads 
to the relation 

2p n 2p K (l- V ) +anP * + 2 M {b) 
This equation of energy conservation may be thought of as a relationship between p t and r\. 
Accordingly, rearranging the terms, we get 

2 _ f( m l- m l)~V [ml + 2pla(l - 77)] 



Pt=V) — x f (7) 

L (l + a)-ar) J 

The minimum and maximum value of r] are obtained by solving Eq. ([7]) for p t = 0: 

2 2 
m^ — m* 

Vrrtin 0] Vniax ~ n \ ~ a (^) 

ml + 2pla 

It is interesting to note that Eq. [8] implies a maximum energy for the neutrino arising from 
the pions of arbitrarily high energy for a given value of a. This maximum energy is given 
by 



E u , max = — £ « « 2.5 x 10 4 GW for a = 1(T 10 (9) 

This limit is noticeable in Fig. [3] 

We next consider the modification of the pion decay matrix element due to VLI effects 
in the Coleman-Glashow model. We begin by writing the pion-decay matrix element in the 
standard form: 

= {t*(Ai)7«(l - 7>K)} UP: (10) 

The symbols fi and in the brackets next to the wave functions are introduced to indicate 
the particles they refer to. The essential change with respect to the standard calculation is 
that we now have 

v{y^)v{y^) =p-m u + a^-p (11) 
After some simplification we have 

\^ 2 =(f^) 2 T (12) 



T vtl = [2{p n ■ p M ) {p^ ■ p v - a(p w -p v )}- ml {p u ■ p^ - a{p u ■ p^)}) (13) 



The following simplifications are relevant: 

■ Vn = g ( m l + m l) ~ a (v 2 pl + Pt) 
{(Pn ■ Pu) - a(p n ■ p v )} = - (ml -ml) +a {p 2 t - T](l - r])pl) (14) 

{{Pu ■ Pp) - a{Pu ■ P^)} = 2 ( m l - m l) - <*W>1 

Accordingly, the decay rate of the pion maybe written as 

r f ( 27r ) 4 |/i#|2x4^ ~ d % d3 P» /-,o 

r -=y ^ |m|5( ^-^-^^^-(2^2^ (i5) 

The integration d 3 /^ is accomplished with the choice of 

PiiL = (1 - ^)Ptt, p M t = -Pw and = -r]p n , dictated by the 6 3 (p n - p^ - p v ) part of the 
integral. Suppressing the constants, we are now left with the integral 

r 7rA1 = / -^—5(E 1T -E il -E u )p 1T d7]dipp t dpt (16) 

where we have written explicitly d 3 p u = p n drj dip pt dpt 

Noting that the Jacobian d ( E -*- E »- E ") = = — an( j that ^ integrates to 27r, 

the decay width of the pion is proportional to 

r„ = r : ^ . — rr (17) 



1 + a) {ml + pi (1 - r]) 2 + p 2 t (r]) } 2 



Note that in Eq. ( 17 ) p t is a function of 77 as given in Eq. ([7]) , and the limits of the integration 



are given in Eq. (J8J). As r] max decreases with increasing p n , the decay probability decreases. 
The other LIV effects are contained in the trace T and the denominator of the integrand in 



Eq. (17). The electronic mode of the pion decay is assumed to be without any LIV effects, 
and contributes about 10~ 4 of the muonic mode for a = 0. r e « 1.2 x 10 -4 r 7r _ At (a = 0) 
and the full decay width of the pion may be written as 

+ r ffe (is) 

Noting that the pion lifetime r is inversely proportional to T, we show in Fig. [T] the 
factor by which the pion life-time is prolonged when a is different from zero at various pion 
energies. 
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FIG. 1. We show here the elongation of the decay lifetime of the pion due to superluminal motion 
of the neutrino, for various values of the parameter a as a function of the pion momentum. All the 
curves are normalized to a = 0, for which r is taken to be ~ 2.2 x 1CP 8 s. The additional factor, 
Ex/kik, is included in the propagation equations, so that for a = we get the standard results. 

B. Simple model for Earth's atmosphere. 

We assume an isothermal atmosphere with the density falling off exponentially with height, 
h: 

p{h) = ^ (19) 

Here we take xo = 1030 g cm -2 and the scale height ho = 7 x 10 5 cm to express the density 
p in g cm -3 . Let t be the length from the surface of the earth along the path of a cosmic 
ray incident at a zenith angle 9, as shown in Fig. [2j The column density x(t, 9) that this 
cosmic ray particle arriving from infinity has to penetrate to reach this point is given by 

x(t, 9) = x sec# e - t/{hosecB) (20) 

so that 

\dt/dx\ = h sec9/x = H/x (21) 
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FIG. 2. The left panel indicates a cosmic ray nucleon incident at a zenith angle 9 that suffers an 
inelastic interaction with a nucleus in the atmosphere leading to multi-particle production. The 
surviving energetic nucleon in the atmosphere and the high energy pion travel essentially in the 
same direction, as do the muons and neutrinos arising from the decay of pions. The path length 
'i' is measured from the earth's surface at the zenith angle 9. The right panel displays the overall 
geometry: the neutrinos and muons enter the detector D at angle 6' ~ 9 for the depth dCR and 
9 > it/2. The zenith angle 9 for downward and upward moving particles through the detectors are 
the same. 

Such an assumption of plane-parallel atmosphere is an adequate approximation up to zenith 
angles of ~ 85°. (For a cosmic ray particle arriving precisely horizontally with 9 = tt/2, the 
maximum column density saturates at w 35 x ). 

C. Propagation of cosmic rays in the Earth's atmosphere. 

At the high energies relevant to the present context, the cosmic rays are incident isotropi- 
cally on the top of the Earth's atmosphere. They interact repeatedly as they descend into 
the atmosphere losing energy through the production of secondary particles, mostly pions. 
These pions are produced with low transverse momenta, ~ 0.5 GeV/c and consequently 
the energetic pions and the leading nucleon propagate essentially in the same direction of 
the primary nucleons. The same is true of the neutrino arising from pion decay where the 
transverse momentum imparted to the decay products has a maximum of ~ 45 MeV/c. 
Accordingly, we assume all the products of the interactions or decay, including muons gen- 
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erated through the interactions of neutrinos underground preserve the zenith angle of the 
parent particle. To proceed, we assume a simple power law for the high-energy spectrum of 
cosmic ray nuclei incident on top of the atmosphere 



f n (E, x = 0, 6) = A_( cm 2 -s-sr- GeV)- 1 (22) 

Here A n is a constant, 7 ~ 1.7, and E is the energy of the cosmic ray particle in GeV per 
nucleon. The nucleons interact inelastically with the air nuclei with an effective mean free 
path, A n of ~ 80 g cm -2 generating pions and other particles. The leading nucleon emerges 
from such collisions with a significant fraction, r] n , of the initial energy. Because of this, the 
reduction in flux of primary nuclei in the earth's atmosphere occurs with a mean free path 
A that is significantly larger than X n . Specifically it can be shown that 

1- < rjl > 

where 



A = n 7 - « 120g/cm 2 (23) 



< VZ >= [\v' n yPWnW n (24) 
Jo 

with P(rj' n ) the probability that the leading nucleon emerges with a fraction 7/ of the initial 
energy. We will encounter similar averages; but we will just use the appropriate averages 
without explicitly showing the angular brackets indicating the average. The flux of the 
nucleons at a depth x in the atmosphere is then; 

/n(^x,0) = ^e- x/A (25) 

These nucleons interact inelastically with air nuclei and generate pions which carry an effec- 
tive fraction 7] n of the primary energy and have an effective multiplicity n n . Thus the rate 
of production of pions of energy E in the atmosphere at column density depth x is given by 

5 ,(£, X ,9) = ^-e-M = -A r .e-VA (26) 

The pions interact and decay in the atmosphere and their spectral intensity is controlled by 
the equation 

df n (E,x) Jl m n \dt(6)/dx\ \ 

- q n {E,x,e) - U{E,x) < — + — } (27) 



dx ' ' 1 A,,- Ecr(a, E) 

Here A^ ~ 120g/cm 2 is the interaction mean free path of the pions, and r(a, E) is the pion 
lifetime elongated due to the LIV (a ^ 0) effects at energy E (see Fig. [T]) and the factor 
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\dt(6)/dx\ = hQsec 8/x is the conversion factor from grammage x to path length at zenith 
angle 6. To facilitate a parallel calculation for the a = case, the factor m n /E is shown 
separately in the decay probability. Defining 

e n = h sec 9 m^/cr^a, E) = Hm n /cT(a, E) (28) 



Eq. (27) may be written as 



-q 7r (E,x,9)-f 7r (E,x)<— + > (29) 



dx " I Xjr Ex 



For sec 6 = 1 and a = 0, ~ 125 GeV. The solution to Eq.( 27) simplifies considerably 
for A fa Avr fa 120 g/cm 2 to yield 



f n (E,x) = A, X e -/A^-(7 + D (^j^) 



(30) 



where A n is a constant. The neutrinos arise through the decay of pions in the Earth's 
atmosphere and the calculation of this production rate involves some subtle considerations. 
Neutrinos of energy E are produced in the decay of pions with higher energy E n = E/rj. 
Yet the effective average value of rj is in itself a function of E n and a. In order to address 
this issue, we show the weighted average, < rj u > = r\ v fa 0.75 r\ max as a function of the 
energy of the neutrino, E v in Fig. [3j For any given values of the neutrino energy, E and 
LIV parameter a, we then can read off the effective mean value of r], and find the typical 
energy of the pion E n = Ejr\ v that generated the neutrino. It is at this energy that we 
should evaluate the pion lifetime, r(a,i? 7r = E u /rj u ), which we write as r v . Note that for 
a ^ there are two values of < rj > that occur; because of the steepness of the cosmic- 
ray spectrum, it is only the larger value of < rj > that contributes significantly to spectral 
intensity of the neutrinos. Defining e v = ho sec 9 m n /(c t v ), and using a similar reasoning 



to that used above in deriving equation (29)and (30), the rate of neutrino production may 
be written as 

q u (E,x) = A,x e-^E-^Sl \ } l - -^f- (31) 



or 



q v (E,x) = A n e^E-^W 1 <f- 



e v + E 



= Q U {E) e~ x/K (32) 
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FIG. 3. The value of < rj u > « 0.75 

Vmax is shown as a function of E u for various values of ex. 
A neutrino of energy arises in the decay of a pion of energy E,,- = E u /rj v . Note that for a^0, 
there are two values of < rj v > for each E u . As the cosmic-ray spectrum is steep, only the larger 
value of r) v contributes significantly. Note that for a given value of a, there is a maximum value 
for the neutrino energy, E u , as stated in Eq. [9} 

where we factored the energy dependent and x dependent terms as Q V (E) and e~^l Xlx . 

The final step in the calculation of the neutrino flux due to the source function q v is to 
include the Cohen-Glashow process of energy loss for the neutrinos [65]. Noting that in a 
single emission of an electron-positron pair through this process, the neutrino loses more 
than 70% of its energy, we treat this process as a decay with an effective lifetime tq given 
by 



r G = \E/c(dE/dx) CG \ = El {cuG 2 F a 3 E 6 } = (33) 

Using the constant ui given by Cohen and Glashow [65], we find t G g ~ 6.5 x 10 _11 s. Keeping 
in mind that in writing the differential equation for the evolution of f u (E,x) we need to 



introduce the factor \dt/dx\ for converting grammage to length, as we did before in Eq.( 27) 
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while describing the decay of pions, we write 



d/„ , ho sec 9 

Qu{E,x) - /„ 



dx xcrc 

= Q^e"^ (34) 

xcr G 



To solve this we let 



9(E,x) = f u (E,x)exp 



x'ctg 



-dx 



n 

X \ CT G 



f u (E,x) [ ) (35) 



where x max = x sec9 and get 



IT 



f = QAE)(^Ve-^ (36) 



This leads to 



(1 \ CT G Fornax H 
/ x^e~ x/K dx 
X m ax J J Q 



Q V {E)K [ / u^Ge~ u du (37) 



Xmax / JO 
\XmaxJ \CT G 

The final step of writing the integral as a complete gamma function follows by noting that 

Xmax 3* ^7r- 



f v (E,x) = g(E,x) 



H 

Xmax \ ctq 



X 

h,Q sec6 



+ 1 (39) 



= ^*"(rj ' B " (,+1, (^j r l^ + V (38) 

At the surface of the earth, x = x maa; and the spectral intensity of the neutrinos is given by 

1 / ~ x ( h sec9 

Detectors for cosmic rays and cosmic neutrinos are placed underground to reduce the back- 
ground due to other particles and gamma rays generated by cosmic rays. Consider such a 
detector, D, placed at a vertical depth, d, as shown in the Fig. [2j The straight line though 
D at a zenith angle 9' « 9 (for 9 < 7f/2) emerges from the Earth's surface at distances L z 
and L n respectively. A theorem in Euclidian geometry yields 

L z -L n = (2R - d)d (40) 

13 



or 

-j2 



L n = (2Rd-d 2 )/L z 
and L z , for 6 < 11/ 2 is given by 

(R - d)cos6 - y/(R - d) 2 cos 2 6 + 2Rd - d 2 (41) 



From this we can calculate L n using Eq. (40). The propagation of the neutrino spectral 



intensities is straightforward if we assume their flux is not significantly depleted due to 
interactions and assume only the Cohen-Glashow process to operate. Accordingly, their flux 
at depth d maybe written as 

f v (E,9,L z ) = f(E,x max ,9)exp[-L z /cT G ] 

f u {E, 6, L n ) = f{E,x max , 6)exp [-L n /cr G ] (42) 

We note that the ratio of these two spectral intensities R Z}n (E) is given by 

R z>n (E) = exp [—(L z - L n )/cr G (a, E)\ (43) 

It is interesting to add a comment here that for most neutrino telescopes operating under- 
ground, the geometrical collecting factor isessentially independent of the hemisphere from 
which the particle arrives, i.e it is the same for downward and upward moving particles. To 
the extent we can neglect neutrino oscillation effects, R z , n (E) will be a good probe of the 
the Cohen-Glashow process. Even though the zenith angle 6, in the Earth's atmosphere, of 
the particles entering the detector after traversing the distances L z and L n are the same, 
there could be a few percent differences in the scale height of the atmosphere at the an- 
tipodal points. When the observational data are averaged over a year, the differences will 
be reduced further. More importantly, the cross section for the interaction of neutrinos 
of energy greater than ~ 10 5 GeV is ~ 10~ 34 (.E/10 5 GeV) 2cm 2 , so that the interaction 
probability across the diameter of the Earth is about ~ 30% and increases with increasing 
neutrino energy as E*. Allowance for this has to be made during the analysis of the data 
while searching for the Cohen-Glashow effect. Alternatively, with stringent bounds on a 
obtained from other observations, the asymmetry in the downward and upward intensities 
may be used to estimate neutrino cross-sections at high energies. The vacuum oscillation 
length in meters is ~ 2.5 E v (GeV) /Am 2 (eV 2 ) so that at the energies of neutrinos under 
consideration this process may be neglected. 
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D. Calculation of the spectral intensities of muons. 

This calculation follows along similar lines as that for the neutrinos and becomes indeed 
simpler when we neglect effects of energy loss due to bremsstralung and ionization in the 
Earth's atmosphere in the region of interest. The muons are generated at a rate 



Here 77^ is the effective average of fraction of energy that the muon receives in the decay of a 
pion of energy E/rj^. As noted earlier in the context of calculating neutrino fluxes, because 
the steepness of the cosmic-ray spectrum, the effective average is ~ 0.75 times the maximum 
fraction. Also noting that the minimum fractional energy carried by the neutrino r] min 0, 
i]^ « 0.75, and is nearly a constant independent of the energy of the pion and the value of 
a. The critical energy is £ w = Hm^/cr^. 

The mean lifetime of the muon r M is ~ 2.2 x 10~ 6 s so that even at ~ 1 GeV its decay 
length is ~ 6 km, roughly equal to the scale height of the atmosphere. Thus in the calculation 
of the spectral intensities of the muons in the atmosphere at energies greater than about 
10 sec 9 GeV, we may safely neglect the decay of the muon. We may also neglect the energy 
losses due to ionization at E > 30 sec 9 GeV. Thus the muon intensity, f^(E,x max ,9), at 
the surface of the earth is given by the integral of the source function. 



III. COMPARISON OF THE THEORETICAL SPECTRAL INTENSITIES WITH 
COSMIC RAY OBSERVATIONS 

The calculations presented in the preceding section indicate that the spectral intensities of 
neutrinos and muons generated in the Earth's atmosphere through pion decay are sensitive 
to the posited level of LIV in the maximal attainable velocities of v^. We illustrate the main 
effects of the superluminal motion of on the penetrating components of cosmic rays in a 
sequence of figures: Figs. [4] - [10] The effect of increasing a on the spectrum of neutrinos is 




(44) 




(45) 
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FIG. 4. The spectrum of neutrinos for various values of a and a fixed value of sec 9 = 10 is 
displayed. Notice that the spectrum steepens sharply at progressively lower energies for increasing 
values of a. For the extremely small value of a = 10~ 20 there is no perceptible steepening even up 
to 10 7 GeV. 



shown in Fig. [4| where the value of sec 9 is fixed at 10, and a is varied in the interval 1CT 20 
to 10~ 6 . Here we can see the neutrino spectra steepening from progressively lower energies 
with increasing a. For the smallest value of a = 10~ 20 , there is no perceptible steepening 
even up to ~ 10 7 GeV. In Fig. [5j we show the dependence of the neutrino spectra on sec 9, 
for a fixed value of a. The well-known enhancement of the intensities at high energies with 
sec 9, due to the increased fraction of pions which interact rather than decay, is reproduced 
even when the superluminal effects are included. The propagation of neutrinos through the 
Earth is exclusively determined by the Cohen-Glashow process, to the extent the neutrino 
interactions with the material of the Earth may be neglected or accounted for. The ratio 



of the neutrino spectral intensities R z (E,9,a) at 9 and 9 + tt calculated in Eq. (43) are 
displayed in Fig. [6j The 9 dependence of the ratio spectral intensities is displayed for 
a = 10" 11 for a set of neutrino energies E = 1 x 10 4 to 1 x 10 5 GeV in the top panel and 
ratio for a fixed E = 10 5 GeV for selected values of a is shown as a function of energy in 
the bottom panel of Fig. |6j 

Similarly we show in Fig. [7] the dependence of the muon spectra on sec 9, for a fixed 
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E v (GeV) 

FIG. 5. The dependence of the neutrino spectra on sec 9, for a fixed value of a = 1 x 10 -10 . The 
enhancement in the intensities at high energies with sec 9 is seen, even with the presence of the 
LIV effects. 

value of a, and in Fig. [8j we show how the muon spectra become steep from progressively 
lower energies as we increase the value of a, for a fixed value of sec 9. We recall that 
our calculation neglects ionization and other losses of energy suffered by muons as they 
propagate through the Earth's atmosphere. The effects of such losses will be to flatten the 
spectra of muons at low energies. We now proceed to compare these theoretical estimates 
with the available data and derive the bounds on the LIV parameter, a. The first cosmic-ray 
observations of muon neutrinos date back to half a century or more, and the instruments 
have progressively increased in collecting power to achieve good sensitivities that we can 
observe cosmic ray neutrinos even up to ~10 6 GeV. Since the LIV effects in the Coleman- 
Glashow model increase with increasing energy, these observations probe very sensitively 
the effects of such violations. 

A. Comparison with cosmic ray neutrino intensities. 

The measurement of cosmic ray neutrino fluxes started with the pioneering efforts of Reines 
et al. [77] and Achar et al. [75]. Progressively the size and sophistication of the detectors 
have improved so much that today we have the spectrum well measured by the IceCube 
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FIG. 6. The asymmetry in the spectral intensity in the forward/backward direction is displayed 
as a function of the arrival direction of the neutrino-induced muons in detectors placed ~ 1 km 
underground, for various values of the LIV parameter a (lower panel). Note that we have shown 
Rz,n — 1 along the y-axis, in order to clearly bring out the dependence on the parameters, (see Eq. 
43 in the text). 



collaboration up to ~ 3 x 10 5 GeV pU]. The observed spectral intensities of neutrinos as 
reviewed by IceCube [10] is displayed in Fig. |9} In the same figure we superimpose the 
theoretical spectra calculated by us for various values of the LIV parameter a, with sec 9 
= 5 representing the weighted average of the intensities over zenith angles of 90° to 180°. 
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E,a (GeV) 

FIG. 7. This figure illustrates the muon spectra generated by the pions in the atmosphere with 
their decay times elongated by LIV effects; the sec 9 enhancement of the intensities at high energies, 
well known in the cosmic-ray field is reproduced, even when LIV effects are present. 

In fact, the effective average value of sec 9 increases with increasing neutrino energy as a 
consequence of the competition between interaction and decay of the pions. At the highest 
energies the mean value of sec 9 is expected to be higher. For comparison, we show in the 
lower panel of Fig. [9] the theoretically calculated neutrino intensities at sec 9 = 10, along 
with the observational data. We normalize all the theoretical spectra to the observed fluxes 
at neutrino energies ~500 GeV. 

We note that the theoretical spectra for a = 10~ 20 ~ fits the observations well. As 
a increases above 10 -14 , the theoretical curves start falling below the observation at the 
highest energies. For a = 10~ 13 , the theoretical curves fall a factor of ~ 300 below the 
observed intensities at ~ 2 x 10 5 GeV, and by progressively smaller factors at lower energies. 
Thus it is safe to conclude that the value of the LIV parameter is less than ~ 10~ 13 . 



B. Comparison with cosmic ray muon spectra. 



Novoseltsev [79] provides an extensive compilation of the spectral intensities of cosmic ray 
muons at energies beyond 1 TeV, and we compare the theoretical spectrum given in Eq. 



(45) with data in Fig. 10 First, we note that for the very small values of a ~ the 
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E M (GeV) 

FIG. 8. The theoretically calculated muon spectra for various values of a, with a fixed zenith 
angle, sec 8 = 10, are displayed. The elongation of pion lifetime due to LIV effects make them 
preferentially interact in the atmosphere rather than decay. The progressive reduction of the high 
energy flux of muons with increasing a is seen clearly. 

theoretical curve reproduces the observations very well, up to ~ 2 x 10 5 GeV. The two 
points at ~ 4 — 5 x 10 5 GeV with large error bars lie above the theoretical predictions. For 
values of a > 1CT 13 , the theoretical curves peel off downwards from the data at progressively 
lower energies. We may thus conclude that these data also constrain the LIV parameter, a 
to be less than ~ 10~ 13 . 

IV. DISCUSSION 

In the preceding sections our aim was to provide a transparent description of the propagation 
of high energy cosmic rays in the atmosphere, resulting in simple analytical expressions for 
the spectral intensities of the neutrinos and muons. Despite the many simplifying assump- 
tions made, our analytical expressions reproduce with adequate accuracy the observations 
and the well-known sec 9 enhancement of the fluxes of the penetrating component at high 
energies. These analytical expressions incorporate the effect of novel LIV physics into the 
propagation of cosmic rays through the Earth's atmosphere and the propagation of neutrinos 
arising from pion decay through the Earth up to the detectors deployed deep underground. 
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FIG. 9. Top panel: The theoretical spectral intensities of the muon neutrinos for various values of 
a are superimposed on the observed spectra for sec 6 = 5, which is taken to represent the weighted 
average of the neutrino spectra at various zenith angles. Note that for very small a = 10~ 20 (~ 0) 
the theoretical fluxes are in good agreement with the observations. Significant deviations appear 
for a ~ 1CP 13 at high energies. Bottom panel: The value of sec 9 has been changed to 10. 



By tracking how the spectral and angular dependence of the muons and neutrinos are 
thereby modified we have placed a strong limit on the LIV parameter: a < 1CT 13 . Several 
remarks are now in order: 



21 




Mars ,1977 
LVD, 12 

Nottingham, 1984 
Moscow University, 1994 
Frejus, 1994 

ufrom pi, K+ prompt, 1994 
ufrom pi,K+VHE 
Artyomovsk, 1985 
Dugapur, 1972 
Baskan, 1992 

1 - Pi, K-muons 

2 - Pi, K-muons + PM(QGSM) 

3 - Pi, K-muons + PM(RQPM) 

4 - Pi, K-muons + PM(VFGS) 
a=1e-20 
a=1e-14 
a=1e-13 
a=1e-12 
a=1e-10 
a=1e-08 



Muon momentum (GeV/c) 



FIG. 10. Differential muon energy spectra from various experiments as reviewed by the IceCube 
collaboration [79] is compared with the theoretical expectation for different values of a and a fixed 
value of sec 8 = 1. The theoretical spectra are normalized at ~ 1000 GeV. Note that for a = 10~ 13 
the theoretical spectrum steepens significantly below the the observed intensities at ~ 10 5 GeV. 
The theoretical spectra shown here are for a fixed sec 8 = 5. 



(i) The electronic decay mode. 

The electronic (e+i> e ) decay mode of the pion is suppressed in the Lorentz Invariant standard 
model by a "helicity" factor of 10 -4 . Furthermore unlike the muons, the energetic electrons 
generated by the cosmic ray u e in the Earth cannot penetrate very far from their production 
site. Thus our analysis and the bounds obtained with an exclusive focus on the muon sector 
do not depend on wether we include or ignore the electronic decay mode. 

(ii) Neutrino mixing. 

It has been pointed out [5U] that the observed oscillations between neutrinos of different 
flavors are strongly suppressed if the LIV parameters a(z/ M ) and a(u e ) differ by more than 
/\m\ 2 jE{y) 2 ~ 1(T 18 for E{u) ~ O(10 MeV) and Am? 2 ~ 1(T 4 eV 2 . Assuming however 
an energy independent common - LIV a{vi) < 1CT 10 seems quite consistent with all other 
data in particular the supernova 1987 upper bound of ~ 10~ 8 on a{v e ) [391 HO]- 

(iii) Charm production. 

A 10 6 GeV neutrino in cosmic rays would arise from the decay of pions of energy ~4x 10 6 
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GeV, which in itself will be generated in the interactions of nucleons of ~ 4 x 10 7 GeV. 
Such an energy for the nucleons correspond to beam energies of ~ 5 TeV in a collider. 
Accordingly, we may expect good data for the production cross section for the very short 
lived charmed and other mesons which may decay readily giving high energy neutrinos. We 
may expect that such processes may start to dominate the neutrino fluxes at these energies 
[8TH83] , making the spectra themselves not so very good probes of the a parameter. On the 
other hand these neutrinos will suffer energy losses through the Cohen-Glashow process and 



the asymmetry parameter R ZyU (E) given in equation Eq. 43 will provide a useful signature 
of the superluminal neutrinos even in this energy region, (iv) LIV for charged leptons. 
Strict SU(2)i gauge invariance in the Standard Model suggests that any non-vanishing 
ct(vi) parameter be associated with an equal LIV parameter for the corresponding lepton. 
Specifically we will then have: 

E(l t ) = [m(l i ) 2 + (l + a(u l )) 2 p 2 } 1 / 2 (46) 

This would imply that asymptotically Ei ~ (l + at)pi so that the muon also becomes equally 
superluminal at high energy. This, in turn, makes prolongation of the lifetime of the pions 
even much more dramatic - since at energies greater than m^/a 1 ^ 2 the pion becomes stable 
and neutrinos or muons of this energy should simply not be produced at all by high energy 
cosmic rays. Note however that if we extend this to a universality of all the a(i) including 
those of the electrons, then the Cohen-Glashow process is kinematically forbidden. An order 
of magnitude estimate indicates that 10~ 10 LIV alpha is (marginally) consistent with the 
precise measurements and calculations of the g-2 of the electron and certainly is allowed for 
the muon and the tau leptons. 

When the Cohen-Glashow process is suppressed, then the leading process would be 

Up -> v + 7 (47) 

which will have a lower rate due to W in the loop. However, the severe limits on this process 
placed by Cowsik, Rajalakshmi and Sreekantan [3T] will apply, 
(v) Connection to GRB's and Supernovae 

A direct 1CT 8 upper bound on the superluminality parameter a(i/ e ) is derived from the 
observed difference in arrival time of about a few hours between the neutrino pulses from 
SN1987a and the optical signature which travel together from the LMC to earth for ~ 3x 10 12 
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seconds [39J. Gamma ray bursts (GRB's), detected by satellites at a rate of about one per 
day seem to originate at cosmological distances of ~500 mega-parsecs ~ 10 4 times larger than 
the distance from the 1987a supernova. Several GRB models suggest that protons should 
be accelerated up to ~ 10~ 20 eV in the fireball and the interaction of these energetic protons 
with the ambient material and radiation should generate pions which eventually decay into 
neutrinos. This motivated the IceCube collaboration to search for coincidences between 
GRBs and upward moving muons - namely energetic muons pointing in the same direction 
[HI]. The good angular resolution (~1 degree) of both IceCube and the satellite detectors 
make for small ~ 0.03 probabilities of one random coincidence even within a generous time 
window of an hour. While no such coincidence has been observed to-date, future observations 
of any coincidences would lead to a 10 4 times stronger bound a(z/ M ) < 5 x 10 -12 , on the 
superluminality of muon neutrinos in the TeV range to be compared with that found for 
lOMeV electron neutrino by using SN1987a data. The bounds set by observing neutrinos 
in coincidence with gamma ray bursts would improve by a further ~100 for a coincidence 
observed within ~1 minute. Amusingly, even using cosmological baseline in putative direct 
timing experiments yields bounds similar to those obtained by our analysis of high energy 
atmospheric muons and neutrinos! 
(vi) G.Z.K and Anita UHE neutrinos. 

Since our bounds improve with the observation of higher energy neutrinos, the searches for 
ultra high energy neutrinos with energies > 10 9 GeV are of special interest. The searches 
for such neutrinos were launched a while ago. The ANITA experiment suspended over 
Antarctica from a high altitude balloon detects the Cherenkov radio emission from neutrinos 
that skim the south pole ice [5]. So far only one candidate event has been found [I] and 
in view of a similar expected background, we cannot establish the existence of such ultra 
high energy neutrinos. Should however future flights provide enough statistics and neutrino 
initiated events of such high energies be definitely detected - the bounds will dramatically 
improve. To see this let us assume that the detected neutrinos are indeed of the GZK 
type namely coming from decays of pions produced via interactions of ultra high energy 
cosmic rays's with the background photons. The neutrinos would then have to travel some 
large distance L ~ 100 Mega-parsecs ~ 10 26 cm to arrive here. Demanding that the mean 
free path for energy loss through the Cohen-Glashow process of these neutrinos exceed this 
distance implies E 5 a 3 < 10~ 26 so that a < 10~ 23 !! 
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V. CONCLUSIONS 



In this paper we have developed a transparent analytical model for the propagation of cosmic 
rays in the Earth's atmosphere that explicitly includes the effects of a superluminal motion 
of the muon neutrino on the decay probability and kinematics of pion decay and energy 
losses suffered by the neutrino through the Cohen-Glashow process as it propagates from 
the production site in the atmosphere onwards to the proximity to the detectors placed 
deep underground. The available observational data to date on the cosmic ray generated 
neutrinos and muons place a bound of a < 1CT 13 . We have pointed out how detectors like 
IceCube may search exclusively for signatures of the Cohen-Glashow process by observing 
the forward-backward ratio of high energy neutrinos arriving at the same zenith angles but 
from upper and lower hemispheres. Keeping in mind that the data published by the IceCube 
collaboration is limited to those acquired in 2009 and earlier years, their full data up to the 
present date should be able to extend the spectra to > 10 6 GeV and the bounds on a to 
~ 10~ 14 . The observations of GZK neutrinos with E ~ 10 9 GeV will push a to well below 
~ 10- 23 !! 
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